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>■ I The production of charmonia in ultrarelativistic heavy-ion colhsions is investigated in- 
cluding two sources. These are a primordial contribution coupled with various phases of 
dissociation, and a statistical coalescence of c and c quarks at the hadronization phase 
^ . transition. Within a schematic fireball evolution, SPS data on J/\l/ production can be 
^ ! reasonably well reproduced. Remaining discrepancies in the ratio are discussed. 

O I Predictions for the J/\l/ centrality dependence at RHIC energies are confronted with first 
' data from PHENIX. The pertinent excitation function of the Nj/q,/Ncc ratio exhibits a 
characteristic minimum structure signaling the transition from the standard J/\l/ suppres- 



Ph. sion scenario (SPS) to predominantly statistical production (RHIC). 



1. Introduction 

In the standard picture of J/\l/ production in heavy-ion collisions [ |T], 0, charmonium 
states are (pre-) formed in primordial A^-A^ collisions and subject to subsequent nuclear 
absorption, Quark-Gluon Plasma (QGP) and hadronic dissociation. Hence, the magni- 
tude of J/\l/ suppression reflects on the plasma effect provided that nuclear and hadronic 
suppression are known and/or small. Recently, an alternative mechanism of J/\E' produc- 
tion in relativistic heavy-ion collisions has been suggested [ |^, |^, |^ . The underlying idea 
is that charmonium states are formed by statistical coalescence of c and c quarks at the 
hadronization transition according to thermal weights, with the total cc abundance being 
determined by primordial parton fusion. 

In the following, we present a combined approach of direct and statistical J/^ pro- 
duction [ constructed within a thermal framework that is consistent with basic 
hadro-chemistry and expansion dynamics at both SPS and RHIC, and also reproduces 
electromagnetic observables at the SPS. In particular, we do not invoke any "anomalous" 
open-charm enhancement beyond expectations from A^-A^ collision scaling. 
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2. Two-component model of J/'^ production 

The first ("direct") component of J/\E' yield is attributed to primordial J/\E''s arising 
from early (hard) parton-parton collisions. A first suppression factor, Snuci is due to 
nuclear absorption evaluated in the Glauber picture with a phenomenological (constant) 
cross-section a^,N — 6.4 mb. Subsequent QGP suppression is based on inelastic scattering 
of thermal partons off the c or c quark within the bound state, using a quasifree approxi- 
mation with account for in- medium modified charmonium binding energies [ ^j. Hadronic 
dissolution cross sections are estimated within a S'f/(4)-symmetric effective theory [ ^, 
The pertinent dissociation rates for both QGP and hadron gas (HG) are convoluted over 



a thermal fireball evolution [ |I0|, |rT| for a heavy-ion collision at given impact parameter, 
yielding charmonium suppression factors Sqq and 5^ {.X= J/"^, ip', x), respectively. This 
leads to the number of direct J/^^'s from primordial production (including feeddown), 

A0% = a'J^ABTMb)Snuc [o.65^/J+^ + 0.0325^^+^ + 0.085^;,^^] . (1) 

The second ("statistical") component of J/\l/ yield originates from coalescence of c and 
c quarks at hadronization. The abundances of charmed particles follow from hadronic 
thermal weights at the critical temperature Tc, with the absolute number matched to 
the amount of charm quarks created in primordial A^-A^ collisions. This necessitates the 
introduction of an effective fugacity factor 7c which is determined enforcing exact (local) 
charm conservation. The statistical J/\t' production, including feeddown, is then given 
by 

A/}/* = IcVh [nj/^S'J'' + E iSniX ^ J/^)nxS^] 7^ , (2) 

where the volume at hadronization, Vh, is taken from our fireball evolution. The factor 
TZ<1 approximates the effects of incomplete charm-quark thermalization (calculated in a 
relaxation time approach), inducing a relative reduction of statistical production. 

The total number of J/\Ef's observed per heavy-ion collision follows as A/j/* = A/'j/^ + 
A/}/^. The model involves essentially two parameters: the strong coupling constant fig- 
uring into the QGP dissociation cross section, as well as the c-quark thermalization time. 

3. SPS 

Our results for the centrality dependence of the J/'^ yield in the P6(158 AGeV)-P6 
system are compared to NA50 data [ |l^ in Fig. [l|. Direct production (dashed line) 
always prevails over the thermal component (dot-dashed line). The latter sets in once 
a QGP starts forming, which, coupled with the rather strong QGP suppression, mimics 
the first drop at around Et—4:0 GeV. The overall good description of the NA50 data 
extends beyond Et > 100 GeV (not shown here) when incorporating transverse energy 
fluctuations [ ^ and a careful assessment of the minimum bias analysis [ , cf . ref . [ 0] . 



In addition, the model also reproduces the NA38 results for the S-U system reasonably 
well. We note that the relative smallness of statistical production stems from not invoking 
any "anomalous" open-charm enhancement. 

Within our two-component model, the \E''/\E' ratio can be extracted without further 
assumptions. The results (Fig. ||) for both S-U (dashed) and Pb-Pb (solid) indicate sig- 
nificant discrepancy with NA38/NA50 data. We believe the deficiency in our description 
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Figure 1. J/^f over Drell-Yan ratio at SPS Figure 2. ^'/^ ratio for P6(158 AGeV)-P6 
for P6(158 AGeV)-P6 collisions. and 5(200 AGeV)-U systems at SPS. 



to reside in the hadronic phase, cf. also ref. [ |T5[. Indeed, an artificial increase of the 



hadronic dissociation cross-section by a factor of 5 clearly improves the agreement 



with the data. 
Z)-mesons. 



Further analysis is required by addressing, e.g., in-medium effects on 



4. Predictions for RHIC and excitation function 



At RHIC, due to much more abundant primordial charm-quark production and due 
to a stronger plasma suppression, the statistical J/"^ production prevails over the direct 
component. This is apparent from Fig. |^ where the number of J/^^'s per binary collision is 
plotted as a function of centrality. The full curve represents our prediction of ref. [ |^, with 
a minor adjustment of the input parameters from pp collisions according to the (central 
values of the) recent PHENIX measurements [ |16| which reported cTj^^ (y^=200 GeV) = 
3.8 ± 0.6(stat) ± 1.3(sys) fih and crc^(-\/i=200 GeV) ~ 650 fih. Despite the large error 
bars, it appears that scenarios involving standard J/\I' suppression only, as well as more 
extreme recombination scenarios, are disfavored by the data. 
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We also note that the probabihstic treatment of nuclear absorption, valid in the low en- 
ergy regime, is not justified at RHIC energies. However, evaluations of nuclear absorption 
at collider energies using a more rigorous quantum-field theoretical approach [ |l^ find 
that the suppression effects are quantitatively very close to the ones in the probabilistic 
treatment. 

The transition between the regimes of predominantly direct (SPS) to statistical (RHIC) 
charmonium production can be further scrutinized experimentally at RHIC in terms of 
an excitation function for the ratio Nj/q,/Ncc, cf. Fig. ^. This ratio exhibits a nontrivial 
shallow minimum structure, a direct consequence of the interplay between suppressed and 
thermal production (assuming no anomalies in open-charm production). 

5. Conclusions 

We have shown that a two-component model for charmonium production in heavy-ion 
collisions, including both "direct" J/^'s arising from primordial N-N collisions subject to 
nuclear, QGP and hadronic suppression, as well as "thermal" J/\E''s emerging from coa- 
lescence of c and c quarks at hadronization, accounts well for J/\E' centrality dependencies 
measured at SPS. A potential discrepancy has been identified in the \E'Y\E' ratio, which 
might hint at shortcomings in the evaluation of hadronic \E'' dissociation. Extrapolating 
our approach to higher energies, the centrality dependence of J /"if production at RHIC 
is in line with preliminary PHENIX data. We furthermore emphasized the importance of 
an excitation function to determine a potential transition from (predominantly) direct to 
(mostly) statistical production, as indicated by an almost fiat (but nontrivial) structure 
of the Nj/ifi/Ncc ratio as a function collision energy. 
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